Collective multipole-like signatures of entanglement in symmetric iV-qubit systems 
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A cogent theory of collective multipole-like quantum correlations in symmetric multiqubit states 
is presented by employing SO (3) irreducible spherical tensor representation. An arbitrary bipartite 
division of this system leads to a family of inequalities to detect entanglement involving averages of 
these tensors expressed in terms of the total system angular momentum operator. Implications of 
this theory to the quantum nature of multipole-like correlations of all orders in the Dicke states are 
deduced. A selected set of examples illustrate these collective tests. Such tests detect entanglement 
in macroscopic atomic ensembles, where individual atoms are not accessible. 
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Correlated macroscopic atomic ensembles [if offer 
promising possibilities in low- noise spectroscopy [2j , high 
precision interferometry [3| , and in the implementation of 
quantum information protocols Experimental char- 
acterization of entanglement - which is the key ingre- 
dient in these applications - has attracted considerable 
attention. Main difficulty in analyzing such composite 
systems with large number, N, of particles is the corre- 
sponding exponential size of the Hilbert space. So, major 
effort is focused on exploring inseparability status of spe- 
cial classes 0, Q of quantum states, confined to smaller 
subspaces of the Hilbert space due to symmetry require- 
ments. For example, a macroscopic atomic ensemble of n 
two level atoms is a collective system of N spin- 1 systems 
(qubits), in a 2 N dimensional Hilbert space H = (C 2 ) 8N . 
However, when the dynamics of the atomic ensemble is 
governed by collective operations - which do not address 
the atoms individually - the atoms within the system are 
completely symmetric with respect to interchange. This 
class of permutation symmetric states, labeled by total 
spin j = ^ (maximal value in the addition of n spin-± 
angular momenta), is restricted Q to a at + 1 dimensional 
subspace n Sym = Sym(c 2 ) 8JV (here ' Sym' denotes sym- 
metrization) . The eigen states {\J - at/2, m) ; -jv/2 < M < 
n/2} of total angular momentum operator j = \ zZ" = i 
where <?; denotes Pauli spin operator of the i th atom, 
span the space. A much simpler analysis of inseparabil- 
ity in symmetric atomic ensembles becomes possible in 
this space. As individual atoms are not accessible in the 
macroscopic ensemble, only collective measurements are 
feasible and any test of entanglement requiring individual 
control of atoms cannot be implemented experimentally. 
For example, spin squeezing [7], i.e., reduction of quan- 
tum fluctuations in one of the spin component orthogonal 
to the mean spin direction below the fundamental noise 
limit jv/4 is an important collective signature of entan- 
glement in symmetric n qubit systems, a consequence of 
two-qubit pairwise entanglement Q . Recently 0] , neces- 
sary and sufficient conditions for pairwise entanglement 



have been formulated in terms of negativity of intergroup 
covariance matrix. These are equivalent to the general- 
ized spin squeezing inequalities [l~oj ] (for two qubit en- 
tanglement) involving collective first and second order 
moments of total angular momentum operator. Further, 
genuine three particle entanglement in symmetric multi- 
qubit systems is shown flpj j to obey inequalities involving 
bulk observables up to third order in total angular mo- 
mentum f. 

In this paper, a family of sufficient conditions to de- 
tect entanglement of atoms in a macroscopic ensemble, 
through collective measurements on the system, is de- 
rived. These are formulated in terms of covariance matrix 
condition involving averages of SO (3) irreducible tensor 
operators ?q(n) of rank k = 1,2, ...n, constructed from 
the total angular momentum operator j and are iden- 
tified via an arbitrary bipartite split of the symmetric 
states of n atoms. The physical significance of these op- 
erators are that they express K = 1, dipole-like; k = 2, 
quadrupole-like correlations etc., among the multiqubits. 
The advantages of this procedure are mainly two-fold: 
one, the simplicity in dealing with a large variety of corre- 
lations in multiqubit systems, and two, the entanglement 
conditions expressed in terms of experimentally observ- 
able signatures associated with the correlations among 
irreducible tensor operators. This elegant formalism en- 
larges the scope of the covariance matrix condition be- 
yond those given in [^, [ll[ . The significance of this new 
approach is substantiated through illustrative examples. 

The SO(3) spherical tensor operators t%(n) 
are constructed such that 

0, 

1 1 31 1 their matrix 
elements in the basis {|A72, a/)} are given by 

(N/2, M'\tq(N)\N/2, M) = sJ2K + 1 C(N/2 K N/2- MQM ), 

in terms of the Clebsch Gordan (CG) coeffi- 
cients 
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These irreducible tensors are orthogo- 
nal, Tt(Tj(JV)fJ,'(JV)t) = (N + i)S KtK ,S QtQl , and the set 
{t§( n )i k = 0,1,2,... N, -k < q < K} forms a linearly 
independent basis of operators in the Hilbert space of 
spin j = n/2 states. Thus a useful representation for the 
density operator of symmetric Af-qubits is given in terms 
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of these operators by 



(1) 



It is completely specified by (a/ + i) 2 -i irreducible tensor 
moments, 



to(N) = Tr[p Bym (AT) t*(N)] = (-!)« t1'{N) 



(2) 



with *o(w) = l because Trp = l. An important composi- 
tion law appropriate for examining multipole correlations 
between two constituent symmetric parts of this n parti- 
cle system, characterized by angular momenta ji = N±/2, 
and n = n 2 /2 (with n = Nt + n 2 ), of the ensemble is con- 
structed in terms of direct product of spherical tensors 

W,JV2) = Wi + 1)W + 1) O^C^^^C'V'^)' 

(3) 

where t«$(N u N 2 ) = Tt[p(N u JV 2 )f~(JVi) ® f*/(AT 2 )] = 

(-1)'+'' /(iVi, A/ 2 ); re = 0,1,..., iVi, re' = 0,1,..., A/ 2 ; 

—re < g < k, —At' < q' < k. Equations ([1]) and ([3]) rep- 
resent the w-qubit symmetric system in two equivalent 
ways and thus the tensor parameters appearing therein 
are related as will be shown presently. These relations 
form the central core of the theory presented here. 

By taking trace over n 2 particles from the composite 
system (using the representation (|3])), we obtain the den- 
sity matrix of Ni particles. We find that the tensor pa- 
rameters tJ(JVi) characterizing the Ni subsystem are given 
by, 

iJ(JVi) = Tr(p(JVi)f»(iVi)) =i^(JV!,JV 2 ); re = 0, 1, . . . , Ni. (4) 

Similar consideration for N 2 -subsystem leads to *j'(N 2 ) = 
t°^,(N 1 ,N 2 ). The second set of relations connecting the 
tensor parameters of the bipartite system given by ([3]) 
with those of , are obtained by using the orthogonality 
property of the tensor operators: 



C'>^) = £ K7'kq {Ni ' N2 ' n) ^ {n) 



(5) 



where |re - «'| < k < k + k'; -k < q < k, and 
f 1k ?'1„{Ni,n 2 ,n) are found to be 

qq' ; KQ l ' 



[*] [f] MMCt^'if;,^) 4 | f 

Ik k k 



(6) 



Here '{ }' denotes the Wigner-9j symbol [12|, [14J and 
[a] = V2a + l. Using the properties of the 9j-symbol and 
the CG coefficients [12j for the special values k = o, q = o 
in (0 and ©, we obtain 



t$(N u N 2 ) =V K {N u N 2 )t*(N), 



(7) 



where v in, n„\ - N i' I (Jv i 1 111 - v ! ,; '• USEE 
o,i,2,. . . Ni. By replacing n 2 -» n 2 - iVi, (where N 2 > iVi) in 



both sides of ([7]) and using ([4]), we obtain an important 
equivalent relation 



t^(JVi) =P«(iVi,iV a - Ni)t£(N 2 ). 



(8) 



The product tensor parameters (Ni,n 2 ) of the bi- 
partite system exhibit a similar relationship with the 
corresponding coefficients t""' («,«') of the « + «' subsys- 
tem. This follows by expressing *q(n) of the compos- 
ite system in the RHS of ^ in terms of « + «' sub- 
system parameters *q(k + «') (obtained by substituting 
ATj = k + k', at 2 = N-(re+re') in ©). Alternately, we can also 
relate t*3'(Ni, N 2 ) to tg(re + re') by choosing Ni - k, N 2 = re' 
in ([5]). Comparing the resulting equations for t* K !(N lt N 2 ) 
and (re, re') and u sing explicit expressions for the as- 
sociated 9j-symbols [12J we obtain, after some algebraic 
manipulation, 



tlq,{Ni,N 2 ) = f(Ni, K ) /(N 2 ,re') t™, (re,re'), 
f(N a , K ) = 



(N a + l)(2re + 1)! 



(re + l)(N a + re + 1)1 (Na - re)! 



a = 1. 2. 



Equations (0, j8|) and (|9|) prove to be significant in 
identifying collective signatures of entanglement in sym- 
metric atomic ensembles, obtained through an analysis 
of the bipartite representation ([3]). 

Consider a set of 2 (2k + i) operators = f~(Ni) ® 
7n 2 and Bl = i Nl <g> f*(N 2 ) (where l Ni = f£(Ni) denotes 
the identity operator). Arranging them as a column £ M 
(corresponding row of operators being £ M = (A^,Bf), 
define the 2reth order covariance matrix for the symmetric 
system as 



i{A 5 -,AC-}. 



(10) 



where A£<*' = f w - (£ (k) ) and {A5< K) , a^ k) } =A£< k) ac< k) + 
Ad K) Af;! re) . y <2K) exhibits a (2k + i) x (2k + l) block 



matrix form, v {2k) = 



A {2 ^(N t ) C (2k) (JVi, AT 2 ) 
C < 2 «)t(Ar l!A r 2 ) S< 2k >(AT 2 ) 

The diagonal blocks A ( "^(Ni) = | ({Alj, Ai",}), 



£T , (JV 2 ) 



■ ({ As" AB", }) correspond to multipole cor- 



relations among the intra-group tensors and the off- 
diagonal block cffi(Ni,N a ) = 5 ({AAJ, ab*,}) comprises of 
inter-group multipole correlations. Explicitly, 

C ( * q ?(Ni,N 2 ) = (-1) 9 ' [^^(W!,^)-*^ ^!,^)*^^^!,^)]. 

(ii) 

Here we focus on the (2re + l) x (2k + i) hermitian cross- 
correlation matrix c (2k) and prove the following theorem: 
Theorem : The cross-correlation matrix c (2re) (iVi, N 2 ) of 
a given rank re associated with any partition (iVi, N 2 ) of an 
Af-qubit symmetric system is necessarily positive semidef- 
inite for all separable symmetric A/-qubit bipartite states. 
The sign of the corresponding matrix c <2k) (k,k), associ- 
ated with the 2k atom reduced system with equal parti- 
tion, suffices to determine that of c <2k - ) (n 1 ,n 2 ), irrespec- 
tive of the partitioning. 

Proof: In the product representation ([3]) with an arbi- 
trary partition (a/i, n 2 ), a separable symmetric A/-qubit 
state has the following structure: 

Paop =^ ProP „(Af 1 )«.p„(Ar 2 ), 0<p w <l; J>» = 1 ( 12 ) 
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where p w (N { ) denotes the density matrices of the 
subensemble of Ni qubits and is expressible as {]]), in 

N i K 

terms of f-(w 4 ): p w (N t ) = £ E ^(^0 t*(N u w). 

K =0q=- K 

In a separable symmetric state (|T2t we have, 

C' (^1 > ^0 = Tr (p scp [(f~ (jVj) ® f (JV 2 )]) 

= Ep^tgCJVi.ti))^,^^). (13) 

Without any loss of generality, we may assume that 
Af 2 > JVi. Now, employing ([7]) we express t hj / (N 2 ,w) — 
[v k (N!,n 2 - JVi)] -1 tg(Ni, w), which leads to an interesting 
form for the matrix c (2k) (7Vi, n 2 ) of (fTTj) in a separable 
symmetric state: 

C i ^, ) {N 1 ,N 2 ) = [V K {Nx,N 2 — JVi)] -1 C (2K) (Afi, JVi) (14) 

Consider an hermitian quadratic form 



: X^ C 



t r^ 2 ") 



(N 1 ,N 2 )X K 



(15) 

with jc" e r( 2 "+ 1 ) being an arbitrary real column vec- 
tor, whose spherical components [12| are denoted by 
= (-1) 9 XT q - In a separable symmetric state (fT2|) it is 
readily seen that Q~ op = [V K {Ni,N 2 - JVi)] -1 n" cp > o, with 
n* ep = c^(n 1 ,n 1 )x k = E to p»E 9 x**ij(isri,^)] a - 
E» Pw £<, X , * 'JC^i, w )\ 2 i a positive semidefinite quantity. 
This proves the first part of our theorem. 

The second part of the theorem follows from ([9]) , lead- 
ing to the result 



C (2k) (N u N 2 ) = /(A' 1 , k )/(JV 2:K )C <2k) («,k), 



(16) 



i.e., the covariance matrix c <2k) (JVi, n 2 ) of an arbitrary 
symmetric system is proportional (with an overall pos- 
itive multiplication factor) to that associated with the 
equal partitioning of a 2k qubit reduced system. (Note 
that the covariance matrix c (2K) (JVi, JV 2 ) could be related 
to its reduced system counterpart c {2K) (n 1 ,n 2 ), which is 
then seen to be proportional to c (2k) (ni,ni); m < n 2 . Fur- 
ther, c (2K) (ni,ni) may be related to c (2K) (n,n): n < m etc. 
This can go down all the way upto c (2k) (k,k). Hence, 
the positivity (negativity) of C C2K) (iVi, n 2 ) has its origin 
in the 2K-qubit covariance matrix c {2k) {k,k,), with equal 
partition (k, «).) □ 

Thus, for any arbitrary (pure or mixed) symmetric 
ensemble of qubits, c (2k) (k,k) < o, for various orders 
k = l, 2, . . . ;2« <JV, serves as a sufficient condition of en- 
tanglement and leads to a family of inseparability condi- 
tions associated with the quantum correlations between 
inter-group tensor operators i.e., (aa^ abJ" 1 ). As the 
product tensor parameters G(k,k) and tj$(re, k), <^)(k,k) 
- specifying the covariance matrix c (2k) (k, «) - are given in 
terms of collective tensor moments t% (N), (see (JSJ) , (0) and 
©), the negativity of the covariance matrix c (2k) (k,k) is 



readily expressed in terms of avera ges of symmetrized 
homogeneous 2k" 1 order polynomials [l3j of j, thus lead- 
ing to a family of multipole-like collective signatures of 
entanglement. 

We now show that the spin squeezing inequality [9|, 
is a consequence of the intergroup dipole correlations ( for 
K = i) viz., c (2) (JVi, jv 2 ) < o: In the 3x3 matrix c (2) (JVi, n 2 ) 
we substitute t K ",'(JVi,iV 2 );re, k' = o, l in terms of total sys- 
tem collective parameters t%(N) (see ©) using explicit 
values 12] for Wigner 9j symbols and CG coefficients. 
Following this by a unitary transformation correspond- 
ing to a change from spherical basis [12j e M , n = ±i,o to 
Cartesian basis e*, i = x, y, z leads to 



U C {2) (N 1: N 2 )U 



N 1 i 
■ A[ Z + VH SS 

1 4 TV 



(17) 



(where i denotes the 3x3 identity matrix; 

\j a Jn + J ic)) - (J a )(3n); s a = (J Q ) and 
Hence, c <2> (ati, w 2 ) < o -<=^ 



jVj w 2 



A — 

JV(JV-l) V <«l+2)(N 2 +2) 

uc {2) {N u N 2 ) c/+ < o v+±ss T < which is the known 
spin squeezing inequality [tj, [l(| , deduced here from any 
bipartite division (Ni,n 2 ). (The covariance matrix condi- 
tion of [11] coincides with that given here only for dipole 
correlations.) 

An important application of our theorem concerns 
the Dicke states [HI \f,M)- -f < m < f. These 
states provide an excellent set of physically relevant mul- 
tiatom symmetric states for illustrating other types of 
multipole correlations. The collective tensor moments 
for the Dicke states are t£(JV) = (f ,M|f^(JV/2)|f ,m) - 
[re]C(f K f ;M0M)5 QO - This corresponds to non-zero 
product tensor parameters t*l' g (N lt N 2 ) (obtained from 
([5])) implying that c (2K) (iVi, jv 2 ) matrix is diagonal (see 
the definition (fTTj) ). Except for |f-,± -f ), which are prod- 
uct states, the covariance matrix c <2k) (k,k) of each of the 
Dicke states is negative, for all orders «, showing the 
quantum nature of multipole-like correlations of various 
orders. Recently (l6j . an experimental scheme to recon- 
struct the spin-excitation number distribution of the col- 
lective spin states ( i.e., tomographic reconstruction of 
the diagonal elements of the density matrix in the Dicke 
basis) of macroscopic ensembles containing ~ 10 11 atoms, 
with low mean spin excitations, has been proposed. Im- 
plementation of such schemes would enable experimental 
detection of collective quantum multipole-like correlations 
in macroscopic assembly of entangled atoms. 

As an illustration of our method, we consider symmet- 
ric mixed states of the form 



_ (1 x)_ 

P (JV + 1) ' ' 



i\\ < x < 1, 



(18) 



where x denotes (JV+i) x (N+x) unit matrix; the states |<fc), 
for i = 1,2, 3, are given by |<£i) = |f , - l), \4> 2 ) = \% ,o) and 
|0 3 ) = ^(|f , f ) + |f )). We find the values of x, - using 
the condition C^^iVi,^) < o - as a function of number of 
atoms, for which p 4 of (fT8|) are inseparable. For the state 
p 3 we find that c (2k) (k,k), are all positive for 2k < n and 
the highest order co- variance matrix c (JV) (iv/2, N/2) < o. 
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consequence is that the Dicke states exhibit quantum 
multipole-like correlations of all orders. Moreover, our 
approach is directly applicable to characterize entangle- 
ment in spatially separated bipartite symmetric atomic 
ensembles, for example, two macroscopic gas samples of 
cesium atoms [17 1. 
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4 6 8 10 12 14 16 18 20 

N 

FIG. 1: Threshold value x m i n , evaluated from the multipole insepara- 
bility condition C^ 2k ^{k, k) < 0, with k — 1 to k — 5 (uppermost curve 
to the lower ones respectively) for the states pi (see (a)) and p2 (see 
(b)), as a function of the number of atoms N. (c): £ m i n obtained from 
the negativity of the highest order covariance matrix C^ N '(N/2, N/2) 
for the states pi (' + '), p2 (' D') and ps (' X '). 



These results are presented in graphical form in Fig. [T] 
From Fig. QJa) and QJb) it is clear that dipole (k = 1) 
quantum correlations lead to x min — > 1 for large n values, 
implying that the mixed states p i; i = 1,2 are separable 
throughout the range < x < 1 in this limit. However, 
higher order multipole correlations are more effective in 
revealing that these states are indeed entangled over a 
larger domain of x. The range of inseparability is sen- 
sitive to the difference n s ~ N g of the number of atoms 
in ground and excited states ( which is N - 1 in \4>i) and 
zero in \<j>v)). Further, from Fig. [ljc) we find that high- 
est order quantum correlations lead to x min — » in the 
large n limit, implying that all the three mixed states 
for i = 1,2,3, are entangled in the range < x < 1, when 

N — > 00. 

In conclusion, we have shown here that SO (3) irre- 
ducible tensor representation provides a powerful method 
to investigate inseparability in symmetric multiqubit sys- 
tems. A family of sufficient conditions of inseparability 
to detect multipole-like collective quantum correlations 
derived here may be useful for experimental character- 
ization of entanglement in macroscopic atomic ensem- 
bles. These conditions are generalizations of the ones ob- 
tained in 0, [ll|] . More specifically, instead of the Carte- 
sian tensor product observables of Ref. [ll|], SO (3) irre- 
ducible tensor observables are shown here to character- 
ize any bipartite division (n 1 ,n 2 ) of symmetric Af-qubit 
systems. The two techniques serve different purposes. 
The first considers groups of qubits to investigate the 
intergroup entanglement, whereas the second describes 
them in terms of spherical tensors involving collective 
w-qubit angular momenta. The latter sheds light on en- 
tanglement among dipolc-like, etc. multiqubit correla- 
tions, which may be physically observable fljj in macro- 
scopic ensembles of symmetric atoms. An important 
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